Generalized rotating-wave approximation to biased qubit-oscillator systems 
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We present the generahzed rotating-wave approximation to biased qubit-oscillator systems in- 
cluding the counter-rotating interactions analytically. Concise analytical expressions for all eigen- 
solutions are explicitly given. Our analytical results for the spectrum of the flux qubit coupled to 
superconducting oscillators can account excellently for the experimental observations. The result- 
ing energy levels for the different biased case agree well with the numerical exact ones in a wide 
range of coupling strength. The dynamics of the qubit are also calculated, the validity of the pre- 
sented analytical approximation is shown to be very well. The present theoretical approach can 
be easily implemented, and successfully applied to the present and near future superconducting 
qubit-oscillator experiments with a biased qubit for accessible coupling strengths. 
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I. INTRODUCTION 

The two-level system (or qubit) and the harmonic oscil- 
lator are attractive components of quantum systems due 
to its interesting application in two-level atoms coupled 
to optical or microwave cavities P, [Ij, superconduct- 
ing qubits coupled to superconducting resonators 
In the early work on cavity quantum electrodynamics 
(QED), the qubit-oscillator coupling strengths g achieved 
were much smaller than the cavity transition frequency 
Lu, g/uj ~ 0.001. All experiments can be well described 
by the Jaynes-Cummning model with rotating-wave ap- 
proximations (RWA) [8]. 

In recent circuit QED, where superconducting artifi- 
cial two-level atoms are coupled to on-chip cavities, the 
exploration of quantum physics has greatly evolved in 
the ultrastrong coupling regime, where the atom-cavity 
coupling strength is cornparable to the cavity transition 
frequency, ff/w ~ 0.1 [l,[i,i,[il[23]. It is evident for the 
breakdown of the RWA and the counter-rotating terms 
are expected to take effect. There have been numerous 
of theoretical studies on the qubit and oscillator system 
finding new phenomena in the ulrastrong couplin g regini e 
pi]| - [ia | and deep strong couphng regime g/u > I |13l.[l4|. 
Recent, Irish J^] has put forth a generalized rotating- 
wave approximations (GRWA) associated with the adia- 
batic approximation and the standard RWA, which works 
well in a wide range coupling regime. However, it has 
not been applied to study biased qubit-oscillator systems 
with a static biased qubit. In circuit QED systems, the 
magnetic energy bias manipulation of the flux qubit is 
associated with the circulating current in the qubit loop 
and the applied magnetic flux 0, [13, I20II. Taking the 
qubit 's bias into account, Grifoni et al. [ll| present the 
Van Vleck perturbation (VVP) theory beyond the RWA 
to treat analytically the two-level system coupled to a 
harmonic oscillator. However, it gives unphysical cross- 
ings in the energy levels in weak coupling region for pos- 



itive detuning. 

In this paper we develope the generalized rotating- 
wave approximations to biased qubit-oscillator systems 
in the ulstrastrong coupling regime, so-called the biased 
generalized rotating-wave approximations (BGRWA). It 
is improved to evaluate a renormalized Hamiltonian with 
the mathematical structure of the ordinary RWA in the fi- 
nite bias case, giving the analytical eigenvalues and eigen- 
functions. It is on the equal foot of the GRWA and can 
recover the results of GRWA for an unbiased qubit. Fol- 
lowing the circuit QED experiment observation in the 
ulstrastrong coupling regime, we provide an concise spec- 
trum analytically using fitting parameters in the exper- 
iment. And we discuss the validity of the energy levels 
of the system by comparing with those obtained by the 
VVP method and numerical exact diagonalization. Fur- 
thermore, we also study the qubit dynamics in the ultra- 
strong coupling regime to confirm the effectiveness of the 
BGRWA. 

The paper is outlined as follows. In Sec. II, the analyti- 
cal solution by the BGRWA to the biased qubit-oscillator 
system are evaluated in detail, giving the eigenenergies 
and eigenstates. And the spectrum is shown with our 
analytical expression using the experiment parameters. 
In Sec. HI, we calculate the dynamics of the qubit in the 
finite bias case. Finally, a brief summary is given in Sec. 
IV. 



II. ANALYTICAL SOLUTION 

The Hamiltonian of the qubit can be written as Hq = 
— {sax + Acrz)/2, where (j^ and cr^ are Pauli matrixes; 
A is the tunneling parameter between the upper level 
I -I- z) and the lower level \ — z) in the basis of ; e is the 
magnetic energy bias related to the circulatin g cu rrent in 
the qubit loop and the apphed magnetic flux |l8l.l20l|. In 
weak coupling regime, where the interaction strength g 
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exceeds the cavity and qubit loss rates, the RWA can be 
apphed and the system can be described by the Jaynes- 
Cummings type Hamiltonian for zero bias as {h — 1) 

H = '^<^z + OJa^a + g[a)(7- + aa+)^ (1) 

where and a are the creation and annihilation opera- 
tors for the oscillator. It can be solved analytically in a 
closed form in the basis |n)| + z) and |n-|- 1)| — z), where 
the qubit states | ± z) are eigenstates of CTz and the os- 
cillator states \n) (n = 0,1,2,...) are Fock states. The 
ground state obtained is \^g) = |0)| — z) for weak cou- 
pling. However, experimental advances are pushing into 
the ultrastrong coupling region, where g approaches to 
the qubit or oscillator frequencies, the RWA no longer 
holds 0, [3 [13 • Thus, the qubit-oscillator counter- 
rotating interaction a^cr+ -I- aa- needs to be taken into 
account. 

Under a rotation around the y axis with the angle 7r/2, 
the Hamiltonian of the qubit-oscillator system including 
the counter-rotating terms reads 

H = -^(jx-'^az+u!a^a + g{a) +a)a^. (2) 

Using the unitary transformation U = 
exp [— ^tJz (a — a^)] , we obtain the transformed 
Hamiltonian H ~ HU — Hq + Hi, consisting of 



Ho = Loa'^a- g'^/uj ~ -az, 



Hi = —^{cTx cosh[— (a^ — a)] + iay sinh[— (a^ — a)]} 



(3) 
(4) 



Recently, many groups have been devoted to solve 
qubit-oscillator systems by different transformations [Tl], 
We now extend the generalized approxi- 
mation by Irish [l^ to the biased qubit-oscillator sys- 
tem. The simplicity of the approximation is based on 
its close connection to the standard RWA. Consequence, 
the terms retained in Hi correspond to energy-conserving 
one-excitation terms, just as the standard RWA. When 
cosh (a'' - a)] is expanded as 1 + -^ [^ (a'f - a)]^ + 



it is performed by keeping the terms 



containing the number operator a^a — n with the coeffi- 
cient Goin) 



Go{n) = (n|cosh[^(at -a)]|n) 

J2 I, ,2\ 



exp(-2.g7w^)L„(4g7t^^), 



(5) 



where L„ are the Laguerre polynomials. The other ex- 
citation terms as well as a^^, a^,..., which are accounted 
for multi-photon process, are neglected within this ap- 

^ a) 



proxiniation. Similarly, by expanding sinh[|^ {a) — a)] 



a' 



a 



the 



one-excitation terms are kept as Fi (n) — aFi (n) with 



the coefficient Fi (n) to be determined. Since the terms 
aFi (n) and Fi (n) involve creating and eliminating a 
single photon of the oscillator, it can be evaluated as 



Fi(n) = 



1 



y/n+l 
5A 



1 1 sinh 



2.9 



1 



(6) 



Since the higher-order terms of Hi are discarded, we con- 
struct an effective Hamiltonian H = Hq + Hi , 



Hq — uja^a 
A 



2/ £ 
A 



Hi = - — [Go{n) - ■q\a,j, - i—(jyFi{n){a^ - a), 



(7) 



where the parameter r? is defined as ?7 = Go(0). 

Since the qubit and oscillator are decoupled in Hq and 
its qubit part can be diagonalized by a second unitary 
u V 



transformation S 



with u 



V2 



1 



-, V 



-^y^l + |, and y — \/ ^ A'^rf. The diagonalized Hq 



IS 



Hq = HqS — ujo'a — g^ 



(8) 



where the tunneling parameter is renormalized by 
i y/e^ + A^Ti^. And the Hi is transformed into 

Hi 



S+HiS 

AajGoin) - Tj] ^ 
2v/e2 + A2ry2 " 



A 



Fi{n){a+ - <j_){a^ - a) 



(9) 



To establish a mathematical structure of the standard 
RWA, the Hamiltonian with BGRWA can be evaluated 
approximately in the form 



H 



BGRWA 



^ 2 + A^4Go{n) - rj] 
va — q UJ ^ — CTz 



2y/e^ + A2772 



(10) 



where the effective coupling strength is renormalized 
into Rr ~ ^F{n), depending on the parameters A and 
g. Combining the transformation, the new Hamiltonian 
with the mathematical structure of the ordinary RWA 
contains the counter-rotating interactions, which play a 
important role in the ulstrastrong coupling regime. 

In the zero bias case, the GRWA by Irish [IJ] works for 
arbitrarily large coupling and for a wide range of detun- 
ing values. We extend the derivation to the BGRWA, 
which can be applied to zero- and finite-bias qubit- 
oscillator systems. As it can be seen from ( [TU|) . the 
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effective Hamiltonian including the counter-rotating in- 
teractions is a generalization of the energy-conserving 
terms i?j.(a^o'_+a(7+) in the usual RWA Hamiltonian ([ij. 
Our approximation exhibits the simplicity as well as the 



GRWA and is expected to extend the range of validity to 
ulstrastrong coupling strengths for a finite- biased qubit. 
One can easily diagonalize the Hamiltonian in the basis 

of I + z, n) and | — z, n + 1) 



H 



BGRWA 



I I 1 \ 2 / s^+A^'jjGoCti+I) 

c.(n+l)-.g 



It is straightforward to obtain the eigenvalues 



2 4Ve^ + A2?72 



■[L„(49V^')-in+i(45V^')] 



±{[ 



2e2 + A27^e-23V"'[L„(V/a;2) +L„+i(45Va;2)] g2^2g-VA 



(11) 



w2(n-f 1) 



(12) 



and the corresponding eigenfunctions 



\lpX) = cos-|n)| + z) +sin-|n-h 1)1 - z), (13) 



= sin-|n)| + z)-cos-|n+l)|-z), (14) 



and b 



In the case of 



where cos 6* = , 

^/<52+4i?,2(ri+l) 
2e^+A^r,[L„(4gVc^^)+L„ + i(4gVc^")] _ ^ 

e = 0, the eigenvalues in Eq.([T2]) are reduced to GRWA 
form 12^. 

The ground-state eigenenergy is 



E. 



BGRWA 



g2 _|_ /^2^2 _ g j^jj^ 



(15) 



with the ground state | — z,0). 

Attributing to the advanced experiment in qubit- 
oscillator systems, there are observations of the spectrum 
in the ultrastrong coupling regime [l^ , which was fitting 
numerically by exact diagonalization in the Fock basis 
[Tsj and in the coherent-state basis 19]. In the exper- 
iment of a flux qubit coupled to an superconducting os- 
cillator, the bias parameter e = 2/p($ — $o/2) with Ip 
the persistent current in the qubit loop, $ an externally 
applied magnetic flux, and $o = ?i/2e the flux quan- 
tum. Fig. [T] shows the spectrum of the system using 
Egs.f [T^ and f [T5|) . with fitted parameters of the exper- 
imental results 5/27r = 0.82, w/27r = 8.13, A = 4.25 and 
/p = 510nA. No significant difference can be observed 
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FIG. 1: Spectrum of the qubit coupled to an oscillator ob- 
tained from Eqs.([T2)) and (I15[) using the experiment param- 
eters in Ref fisll. 



between our results and experiment ones jl8| . It demon- 
strates that our BGRWA will have a remarkable useful 
application in the future experiment since large strong 
coupling strength will be accessed. 



To the best of our knowledge, the analytical expres- 
sions Eqs.([ni)-([l5l) for the present BGRWA is the most 
simple one among the existing analytical works. It should 
become potentially effective tools in the study of the su- 
perconducting qubit-oscillators where the biased param- 
eter is mainly manipulated, if the validity is verified in 



4 



4 




numerical 






- - VVP 


2 




• our method 




















-2 








s=0.1 




-4 


■ A=0.5 





-4 



0.0 
4 



0.5 



1.0 



1.5 





numerical 





VVP 
• our method 


















£=1 








■A=0.5 





0.0 



0.5 



1.0 




1.5 



FIG. 2: Energy levels En as a function of coupling strength g for different bias e = 0.1, 0.5 and 1 (from top to bottom) with 
the tunneling parameter A = 0.5, 1 and 1.5. We compare the eigenvalues in Eq. H12p (solid dot) with those obtained by the 
numerical exact diagonalization method (solid line) and VVP eigenvalues in Eq. ([16]) (dashed line). We set a; = 1. 



general. To show the validity of the BGRWA in the bi- cussion on the energy spectrum and consider eigenvalues 
ased qubit-oscillator system, we present a detailed dis- obtained by the VVP method, which is written as [ll| 

I 



(to + -)a; - — + - > ( 



Dl 



k=0,k^m± 



Dl 



^ £ + (to — k)Lj £ + (fc — n)u 



\ 



[e - + A2 ^ ( 



D2 



^nk 



k=0,k^m±l 



£ + (to — k)uj e + {k — n)ui 



)12 +4£)2 



{n = m + l,l>0) 



(16) 



where Z?m„ = ■j(— l)™(2g)"~™e~^^^i/^LJ^^™(4g^/cL;^). The mth eigenvalues is a mixture of the oscillator 
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levels m and /. Note that the ambiguity of the value of 
I, which should be selected by experience to give better 
results. It works well for large bias e and strong cou- 
pling. By comparing with our analytical eigenenergies in 
Eq. (jl2p . our approach can be more easily implemented. 
The detailed comparison with our approached is present 
in the following for different coupling strength g and 
detuning parameter A/w. 

Fig. [2] shows the first eight energy levels against the 
coupling strength g for three different values of bias e and 
tunneling parameter A for energy scale uj = I. For neg- 
ative detuning A = 0.5, our analytical method and the 
VVP show good agreement with the numerical results 
from weak coupling regime to strong coupling regime 
for different bias values e — 0.1, 0.5 and 1, as shown 
in Fig. (^fleft panel) . At resonance A = 1 (middle panel) , 
our analytical solutions agree well with the numerical 
results at g < 0.5, which is just the present accessible 
coupling regime g < 0.12 and will be first possibly 
accessible regime in the future experiments. Just at this 
practically interesting coupling regime, the VVP results 
deviates considerably. The deviations from the numeri- 
cal results in the intermediate coupling regime is due to 
the dominated infiuence of the higher order terms ne- 
glecting in the transformed hamiltonian in Eq. ([7|) where 
more photon excitations play roles. In the case of posi- 
tive detuning A — 1.5, there is an improvement by our 
analytical results in weak coupling regimes with its com- 
parison to the VVP ones, as shown in Fig. [5] (right panel). 
Especially, for e = 0.1 and 0.5, the results of the VVP in 
the weak coupling regime is qualitatively incorrect with a 
unphysical crossing, but that of our BGRWA solutions is 
still in agreement with the numerically exact one. There- 
fore we can say that the BGRWA approach, which takes 
into account the effect of counter-rotating terms, is quite 
simple and can give accurate analytical solution in the 
biased case. 



III. DYNAMICS OF THE QUBIT 

In the original Hamiltonian ( [5]) , the excited wavefunc- 
tions without RWA can be evaluated by unitary trans- 
formation l^*^) = U'^ S^\ip^) in the following 

l^+^BGWA ^ {ucos^\n)^+vsm^\n + l)^)\+x) 

9 9 
+ (vcos -\n) -3 — ttsin-|n -I- 1) -g)| — x), 

^^-^BGRWA ^ (ysin^ln)^ -i;cos^|rH-l)^)|-|-:E) 

6 6 
+ (^sin -|n)^ — ucos-\n -f 1)^)| — x), 

(17) 

where the qubit states | ± x) are eigenstates of <Tx 
and the oscillator states = e^^/'^'""'^^-' |n) are 

displacement-shift Fock states, so-called coherent states. 



And the ground-state wavefunction is obtained by 
\xi;g)BGRWA ^ „|e)e-9/'^(«-«"^)|0) - u|c,)e5/"(°~°') |0) . 

(18) 

We now contrast the dynamical evolution of az{t) to 
present the validity of the analytical solutions. The initial 
state is assumed as \<f{0)) = \ t)|0). With the eigenba- 
sis {\^n)^'^^^^} and eigenenergies {E^'^^^^}, the dy- 
namical wavefunction of the Hamiltonian without RWA 
( can be expressed as 

\ipit)) -e-''«V(0)) 

= ^e-*^^""""'^|vl/„)«^™(*„|^(0)). 

n 

(19) 

We calculate the time evolution of {<Jz{t)) = 
{ip{t)\a z\^{t)) using our method presented in this paper 
without RWA for the biased case. The expectation value 
of {iTz{t)) is plotted in Fig. [3] for e = 0.1 (upper panel) 
and \/0.5 (lower panel) for different coupling strength 
g = 0.1,0.2,0.5. For comparison, the results of exact 
diagonalization and those of VVP are also shown. One 
finds that the analytical dynamical evolution agrees well 
with the numerical exact ones, much better than those 
obtained by VVP. It follows that the contribution of the 
counter-rotating interaction is reasonably considered in 
the BGRWA analytical solution. Thus, the BGRWA can 
be applied validly in the dynamical evolution in a wide 
range of coupling strength. 

IV. CONCLUSION 

We have presented the analytical solution to the qubit- 
oscillator system without the RWA for the finite-biased 
qubit by the BGRWA. The simplified approach takes 
into account the counter-rotating interaction well but still 
keeps the simple mathematical structure of the ordinary 
RWA. It facilitates us to give the analytical eigenfunc- 
tions and eigenenergies simply, which can also recover the 
results of GRWA for zero bias. Our approximation has 
extended the validity to the ultrastrong coupling regime. 
The spectrum obtained by the analytical expression are 
in good consistent with experiment observations. And 
we have compared with energy levels with the VVP and 
numerical results, exhibiting a wide range of validity to 
large coupling strengths g < 0.5. One finds that there 
is an improvement of the energy levels of ground and 
lower-lying excited states by compared with VVP results. 
Especially, our analytical energy levels fit well with nu- 
merical results in the positive detuning region, where the 
VVP is invalid. Moreover, our BGRWA dynamical evolu- 
tions quantitatively agree well with the numerical exact 
ones for weak to ultra-strong coupling. The present ana- 
lytical approach is simple to be implemented and should 
be very useful for the present superconducting qubit- 
oscillator coupling systems with g < 0.12. More over. 
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FIG. 3: Time evolution of az{t) for £ = 0.1 (upper panel) and \/0.5 (low panel) for different coupling strength g = 0.1, 0.2, 0.5 
on resonance. The dashed line is our BGRWA calculation. The solid and the dotted line are the numerically exact result and 
the VVP one, respectively. The evolution starts with a vacuum state |0) and an excited spin |e). 



there is still wide room g < 0.5 in BGRWA to accu- 
rately account for the future experiments in the same 
field where the coupling strength is further increased to 
explore novel phenomena beyond the RWA. Finally, it is 
suggested that the present approach can be extended to 
more complicated problems, such as biased spin-boson 
systems. 
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